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Several dynamical symmetries of the Dirac Hamiltonian are reviewed in a systematic 
manner and the conditions under which such symmetries hold. These include relativistic 
spin and orbital angular momentum symmetries, S*0(4) x SU a {2) symmetry for the Dirac 
i^« Hydrogen atom, SU(3) x SU a (2) symmetry for the relativistic simple harmonic oscillator. 

J> ' The energy spectrum in each case is calculated from group-theoretic considerations. 
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I. INTRODUCTION 

O j 1 The 1/r potential of force provides the underlying dynamics for the Kepler problem in classical 

mechanics, hydrogen atom in quantum mechanics and one gluon exchange potential in perturbative 



> 

o 



X 



quantum electrodynamics (PQCD). This and simple harmonic oscillator are one of the few physical 
systems that we know how to solve. Further most of more complicated systems can be studied by 
representing them as a collection of harmonic oscillators with different frequencies and amplitudes. 
Not only the above two systems are solvable exactly both in classical and quantum mechanics, 



they also provide the realistic models to study dynamical symmetries as distinct from geometrical 



symmetries. Symmetries do play an important role in the progress of physics. Once they are known 
for a physical system, many of the properties of that system can be established in general terms 
without actually solving for the underlying dynamics. 

As stated above symmetries are of two types: (i) Geometric, the well known example of which 
are space-time symmetries, e.g. rotation (ii) Dynamical, where the underlying dynamics manifest 
some symmetry. This is well illustrated both by the hydrogen atom and harmonic oscillator, which 
are rotationally invariant but also have a well known feature that their orbits of motion close on 
themselves in classical mechanics. This shows that in addition to angular momentum (conservation 
which is due to rotational invariance), there are additional constants of motion. For the hydrogen 

n 

atom this constant of motion, as is well known, is provided by the Rung-Lenz's vector [1J. For the 
harmonic oscillator this is provided by quadrupole moment second rank tensors [2]. 
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For the hydrogen atom in Schrodinger theory, energy levels are given by E n = —7r%, where n 

is the principle quantum number, and a is the fine structure constant. For a given I, n > I + 1; 

I = 0, 1,2, ....n. The {21 + 1) degeneracy with respect to magnetic quantum number is due to 

n-i 
rotational symmetry of the potential — — . For one energy value E, there are Yl (2Z + 1) = n 2 

1=0 

different possible eigen-functions; such a degeneracy occurs only for 1/r potential of force. Usually 



a degeneracy is associated with a symmetry and it is known 



l|] that in this case there is an external 



symmetry 



R,H 



0,with H = ^ — j, and R = -^(p X L - L x p) — jr. The operator R given 
above in Quantum Mechanics follows from the Lenz's vector in the classical Kepler problem by the 
correspondence principle. The orbital angular momentum L and K = \ -^R {K is hermitian when 
acting on eigenstates of H with-negative energy eigenvalues- the ones in which we are interested; 
\i~Tfi commutes with R and L) generate 50(4) algebra, which is isomorphic to SUm(2) xSUn{2) 
generated by M = \ , N = —^ — ■ This symmetry leads to n 2 degeneracy [1[. 

Similarly for the non- relativistic harmonic oscillator, the energy spectrum shows degeneracies 
in addition to those which arise due to rotational invariance. As is well known, the energy spectrum 
En = 2^w(2iV+3) depends on the oscillator quantum number ./V = 2n+l, where n > is the radial 
quantum number and I is the orbital angular momentum. Thus all states with I = N, N — 2, ....0 or 
1 have the same energy. These degeneracies are produced by an SU(3) dynamical symmetry |3j]. 

Obviously the above symmetries are broken in the relativistic quantum mechanics, where a spin 
1/2 particle satisfies the Dirac equation. As is well known this is because the spin-orbit coupling 
leads to splitting of the energy levels [see below]. Since spin is purely a quantum mechanical 
concept with no classical analogue, the guidance from the corresponding principle is not there to 
obtain dynamical symmetries for the Dirac Hamiltonian. As such it is highly non trivial to discuss 



under which conditions the larger symmetries are obtained in the Dirac Hamiltonian. This has 

n n 

been answered for harmonic oscillator in [4| and following the method of [4|, for the Hydrogen 



atom in [5|. The purpose of this paper is to review the previous work mentioned above and to 
discuss the symmetries of the Dirac Hamiltonian in a systematic way, using the Dirac algebra and 
for harmonic Oscillator SU(3) algebra generated in the Gell-Mann basis [6( and point out which 
physical system might have such symmetries. In fact we first show that the harmonic oscillator 
provides a simple realization of 5(7(3), which has been used in particle and nuclear physics. The 
energy spectrum in each case has also been calculated from group-theoretic considerations. 



II. DYNAMICAL SPIN AND ORBITAL ANGULAR MOMENTUM SYMMETRY FOR 

THE DIRAC HAMILTONIAN 

The Dirac Hamiltonian is given by 

H = a-p + Vy (r) + /3m (1) 



where for the Hydrogen atom Vy (r) = —j and is the time component of electroma 
j^i ^y^ ^ _ 0,1,2,3). The hydrogen atom in Dirac theory is exactly solvable [7 
levels are given by 



netic potential 
and energy 



„ 2 



^—n + i^xpi"' 2 



where 



5j = {j + \)-[{j + \?-c?\^ (2) 

and j = 0,1/2,1.... The above expression manifests the fine structure and spin-orbital splitting. 
The states with same n but different j e.g. 2pi/2 and 2p 3 / 2 are now split and the splitting is in 
agreement with experiment , a great triumph of the Dirac theory. The states with same n and j 
e.g. 2si/2 and 2p 1 / 2 are still degenerate, the so called Lamb shift which needs radiative quantum 
corrections provided by quantum electrodynamics (QED). The SO(4) symmetry for the hydrogen 
atom in Schrodinger theory is thus broken in relativistic theory. Usually in the non-relativistic 
limit i.e to order (p 2 /m 2 ), the relativistic theory reduces to non-relativistic theory but in this 
case the non-relativistic theory has a higher symmetry, namely, SO (4) which is not maintained 
by the relativistic theory. The question then arises under what circumstances for hydrogen-like 
system (by this we mean systems with 1/r force potential), the SO(4) symmetry is restored for 
the Dirac Hamiltonian, although it will be contrary to experimental observations for the hydrogen 
atom. Then why should even this question be considered? This is because there are other physical 
systems which show such an approximate symmetry as we will discuss shortly. To answer the above 
question, let us introduce a Lorentz scalar potential Vs (r), then the Dirac Hamiltonian becomes 

H = a-p + V v (r) + f3(V s (f) + m) (3) 

and the corresponding Dirac equation is 

[n°(d + iV v {r)) + iy% -m- V,(r)]V = 0. (4) 



where -y l = /3a*,7° = /3 

If we multiply on the left by 

[i~f (d + iV v {r)) + ij j dj + m + V s (r)]. 

Then for stationary states (Jr — > iE) , the Dirac equation becomes 

„o„i 



where 



[V 2 + V? - ^ - 2EV V - 2mV s + 17 V [$, V v ] - i 7 '[$, V s ] + (E z - m 2 )]* = 0. (5) 



dV 
ft, V»] = £• W 



Thus if 



^s(^) = V(l>) + C/ 5 (7) 

where TJy and C/s 1 are constants, then since in the non-relativistic limit 7 — > 1 + 0(-^),7* —> 
1 + 0(£-), E — > e + m, the terms of order ^ cancel out in Eq. (5). As a consequence absorbing the 
constants Uy and Us with re-definition of E and m, the above equation reduces to the Schrodinger 
equation, which as already seen has 50(4) symmetry, if V is spherically symmetric. 

Indeed it has been observed that the Dirac Hamiltonian (3) is invariant under a spin symmetry 



y, [h, 



s 



0, provided that conditions (6) or (7) are satisfied. Here the generators S form the 



spin SU(2) algebra and are given by 



S=\ I (8) 

Up s Up 

where s = a/2 are usual spin generators and u p = — is the helicity unitary operator. It is easy to 
check that 

[Si,Sj] = i €ij k S k (9) 

In the Pauli representation of Dirac matrices, 

P= I , a. = U ) , 7 5 = -ia l a 2 a 3 = (10) 






By introducing £ 



a 

8 



, one can write the Dirac Hamiltonian (3) and the spin operator (8) 



as 



H = -f 5 [3 S • p + VV (f) + p(y s if) + m) 



(11) 



«-k 



/5S + (l-/3)S-pp-i 



(12) 



Then, since [/3, 1 - /3] = 0, /3, £ 



0. 



7 5 ,S 



0, [ 7 5 ,/3] =0and 



£ • p,T, • p p l 



H,S 



lb 5 '/ 3 ] 



Z-p, £' 



+ 4 [7 5 , (I"/?)] 






+ 



(1-/3) 



W (r),(l-/3)^^ 



+ 



/3(l-/3) 



/r 



pi 



I 7 5 /3 2 dV - ^7 5 /3 2 p* + <LJ±V 



(Vv(f)-V s (r)),^4- 

pZ 



(13) 



Thus 



H,S 



if 9*VV (0 = diy s if) or V~y (r) = V> (r) + U i.e. if the conditions (6) or (7) 



are satisfied. Further for spherically symmetric potentia 



Dirac Hamiltonian has an additional invariant algebra 






s, V v if) = V v (r), V s (r) = V s (r), the 



[9],[ 



H,L 



where 



L 



I 



Up I Up 



(14) 



and I = r x p is the orbital angular momentum. One can write (14) as 

Il-'E-pp/p 2 



L = f+<I^> 



Then 



VV (r)J 



0. 



VsirU 



and it follows as above that 



H,L 



0, [Li, Lj\ — i €{jk Lk 



(15) 



(16) 



Thus S and L are separately constants of motion while s and I are not, but s + /is. 

Finally the Dirac Hamiltonian, which has the above relativistic dynamical spin and orbital 
angular momentum symmetries, is 



H = 7 5 £ • p + (1 + /3)V(r) + (3m 



where constant U can be absorbed in the mass term m. We note that 



(17) 



(H + m) = 7 5 S • p + (1 + /?)(V + m) 



(18) 



H 2 -m 2 =f + 2(1 + /3)V(m + V) + 7 5 (1 - P)Vt ■ p + (1 + 0)VE • p 



(19) 



There is another symmetry the so called relativistic pseudo spin and orbital angular momentum 
symmetry [9J, which is obtained from the one considered above by making the 75 transformation 
and m — > —m 



H -> H = 7 5 #7 5 = 7 5 S • p + (1 - /3)V(r) + (3m 



s^ = ^ s ^ = h-pt + (i + p)^} 



Z^-U^V-lj*] 



[H,~3] = 0,[H,t]=0 



(20) 



The first of Eqs. (20) implies Vy = — Vp+ constant for pseudo-spin symmetry. The pseudo-spin 



degeneracies have been observed in nuclei 



lOj. In fact relativistic mean field representations of the 
nuclear potential do have this property |lll | namely Vy — — Vs- 



We now consider the question wether the property Vy( r ) = Vs( r ) + U arises in QCD. In the 
(Qq) or (qQ) bound states meson spectroscopy, where Q is a heavy quark c or b and q is light 
quark u,d, or s, the spin-orbital splitting is seen to be suppressed experimentally (121 ]. In a (Qq) 
system, where Q is considered to be infinitely heavy, the spin Sq of the heavy quark is decoupled, 
it is then natural to combine the angular momentum of light degrees of freedom j = L + S q with 
Sq to give J = j + Sq for the bound (Qq) system. Thus for the p states we get two multiplets, 
one with j = 3/2 and other with j = 1/2 which for Q = c, and q = u or d are [and similar one for 
B meson when Q = b] 




1 



[D*(1-),D(0-)] h , 

[ J D 2 *(2 + ),^ 1 (1+)] J= 3, 
■> 2 

and [D*(l+),A,(0+)] i= i, 
The splitting between j = | and j = \ multiplets is due to spin-orbit coupling L ■ S, 



T^> 



as in 



hydrogen atom) while the hyperfine splitting between the two members of each multiplet arises 



from the Fermi-contact term S q -Sq, the spin-orbit coupling term (Sq+Sq)- L and the tensor terms. 
The splitting [12], for D mesons, between Dl(l+) : 2422.3±0.6MeV and D* 2 {2+) : 2462.8±1.0MeV 
is WMeV; for B mesons, between Bi(l + ) : 5723.4 ± 2.0 and B%(2 + ) : 5743.9 ± 5 is 20MeV; for 
B Sl (l + ) : 5829.4±0.7 and B S2 (2 + ) : 5839.7±0.6 is WMeV. We notice that the spin orbit splittings 
seem to be suppressed. A measure of this suppression [l3j] is the parameter 

_ P3/2 ~ Pl/2 

(4p 3/2 + 2p 1/2 )/6 - s 1/2 

which for the experimental data shown above is of order 0.07 both for D and B mesons. We also 

note that suppression generally increases with the increasing mass of Q. It has been suggested 13] 



that the dynamics necessary for the approximate relativistic spin symmetry discussed above in the 
heavy-light quark system may be possible in QCD. 

III. SO{4) x SU a {2) SYMMETRY FOR A COULOMB LIKE SYMMETRY POTENTIAL 



In this section we discuss the generators of 50(4) symmetry, which the Dirac Hamiltonian (3) 
possess. Since in the non-relativistic limit j3 — > 1 + 0(p 2 /m 2 ) and 7 5 — > 0(p/m), the natural 
generalization of the Lenz's vector in the Schrodinger theory, for the relativistic case is 



where 



where now 



2mR ->• 2mr 



2mr = (! + /?) f(r)f+ A + (1 + /3) >y 5 g(r)rL -p+(l- p) 75S ■ pg(r)r 



A = px L — L x p 



(21) 



(22) 



L is the relativistic orbital angular momentum defined in Eq. (14). Since H involves 7 5 S • p, 
therefore T should involve such a term and the second term in square brackets appears to make 
the operator hermitian. The functions f(r) and g(r) are to be determined from 



H,T 







(23) 



Now L commutes with H and also with V and /3m and therefore it also commutes with 7 s S ■ p , 



it follows that 



7 5 S -p,A 



since p commutes with S • p. Thus using (S • p) 2 = p 2 , 



2m 



H,T 



S • p, f(r) f - 2(V + M)g(r) r 



+ (l + [3)[V,A] + (l + /3)[p 2 ,g(r) 



+ 7 5 /3 



E-p, f(r)f-2(V + M)g(r)f 



(24) 



The condition (23) gives 



and 



f(r) = 2(V + m)g(r) 



(l + /3)[y,A] + (l + /3)[p 2 , 5 (r)r] =0 



Now using Eq. (15) and the fact that (1 + /3)(1 - /3) = 0, 



(25) 



(26) 



(1 + P)[V,A] = (1 + 13) 

= (1 + /9) 



V,px I — I x p 
idV r 



r dr 



r x I — I x r 



(27) 



On the other hand 
[p 2 ,g(r) f 



-2% 



g{r) + p- 
or 



- *#0 r 

r or 



r xl — I x r 






(28) 



It is important to point out that so far we have made no commitment to the form of potentials 
V(r) and g(r). The condition (26) is satisfied if 

g(r) = V(r) 



and 



dV 

V(r) + r— - = 
or 



(29) 



which gives V(r) =constt. - i.e, the Coulomb potential — — . It is the constraint (26) which forces 
the relations (29) and as a result V(r) has to be the Coulomb potential. Thus 

2mf = 2(1 + 0)V(V + m) f + A + F 

where 



^ = 75 



(I- /3)V f E-p + h.c. 



(30) 



In order to find j^ir [P, D 7 '] , we note that 

[A\A->] = -4e ijk p 2 L k 



(31) 



[(1 + P) 2 V{V + m)x\ A j ] -i++j = -4(1 + {3)V{V + 2M)ie ijk L k 



(32) 



where we have used, r^f- = —V(r) for V = — — and that [c.f. Eq. ([15 



Further 



(1 + f3)l k = (1 + (3)L k 



[A\ F]-iHj = -Mi 3k 7 5 (1 - P)V E ■ p + (1 + 0) E • pV 



[(1 + /9)V(V + m)x\F 3 ] -mj = 



L* 



(33) 



(34) 



[F\ F 3 ] = -4(1 + pyi k L k V 2 

Collecting the various terms and using Eq. I|19]) with V = —a/rwe see that 

tf 2 -m 2 



[r',r>] 



_ ie ijk L k 



m" 



Defining K* 



™_ H -iX % ■, we have finally, 



[iT,^] = ie ijk K k 
[L\L j ] = ie ijk L k 



(35) 



(36) 



(37) 



which generate 50(4) algebra. Further S l commutes with L 3 as well as with H and K 3 . Thus 
defining 



M l 

N 1 



2 



we see that 



[M\M 3 ] = ie ijk M k 
[N\N 3 ] = ie ijk K k 
[M\N 3 ] = 



(38) 



Thus the invariance group for the Dirac Hamiltonian (13) for the Coulomb potential of the hydrogen 
atom is SUm <8> SUn (8) SU a (2), where SU a (2) is the group generated by S l given in Eq. (12). 
The energy spectrum can now be easily determined 

2 (H + m) 2 



o H 2 - m 2 , o „>. 

T 2 = = (L 2 + 1) + 

mr 



a 



m'- 



(39) 



10 



and r • L = L ■ F = 0, implying K ■ L = so that M 2 = TV 2 . In terms of Casimir operator, M 2 , we 
can write Eq. (|39p as 

4(m 2 -H 2 )M 2 = H 2 -m 2 + a 2 (H + m) 2 (40) 

Now since M obey angular momentum commutation relations, M 2 has eigenvalues m(m + l), where 

m can take on the values 0, 1/2, 1, • • • . It is customary to use j for m, then Eq. (|4"0~]) . gives the 

energy eigenvalues 

,_, 4n 2 - a 2 . , , , 

^ = -4^W (41) 

where n = 2j ' + 1, i.e. the energy spectrum is determined only by the principal quantum number 

n and states in the different j values are degenerate showing no spin-orbit splitting. In fact Eq. 

(41) gives for e = E — m 

e n = -m— + 0(a ) 
2n 

which agrees with the energy spectrum for hydrogen atom in Schrodinger theory. 

The Dirac equation with vector and/or scalar Coulomb like potentials 

ay 



V v (r) 



r 



s 



Vs(r) = -^ (42) 



is exactly solvable [8(] and energy spectrum is given by 

E = 2+7 xv {~ ava s ± ( n ~ S ^ i a v - 4 + («- ^) 2 ] 1/2 } (43) 



where 

1/ 



Sj = u+h 



(j + \? - (* 2 v - 4) 



(44) 



For ay = as, this reduces to Eq. (J4T1) . 



IV. SU(3) SYMMETRY FOR THE RELATIVISTIC HARMONIC OSCILLATOR 

" The career of a young theoretical physicist consists of treating the harmonic oscillator in ever- 
increasing levels of abstraction." - Sidney Coleman 

In this section we discuss how harmonic oscillator provides a simple realistic model to introduce 
SU(3) symmetry [2| in a language which is more familiar. 



11 



In non-relativistic quantum mechanics, the harmonic oscillator Hamiltonian 



tt ± r 2 . 2 2 2i 

H = - — \p +m u r \, 
2m 



which is symmetric in p i — > x, commutes with the quadrupole moment operator Q 



'J 



Q 



ij 



^X%iS\ , ( JL^ _!«?'„ 



m l u l [ x % x 3 - -b %:s x l ) + ( p l p 



;S»jf 



(45) 

[i,j = 1,2,3] 

(46) 



Note that, being symmetric in i and j as well as traceless, Q 13 has five independent components 
while the orbital angular momentum I 13 = x % p> — x J p l , being antisymmetric in i and j, has three 
independent components. In order to go from tensor basis to the Gell-mann basis, we introduce 



F n 



2 2— T^a— , l l—T^a— —T^a — \ , _tA _ 

m uj x — x H — muj x — p — p — x + p — p 
2 2 V 2 F ^2/ F 2^ 



(47) 



where a = 1, ....8, and x and p are column matrices [belonging to representation 3 of SU(3)] 

fx x \ ( p l \ 



,P 



\ X 3 j 



(48) 



\P 3 / 



The superscript T denotes transpose so that x and p are row matrices. The tensor and Gell-Mann 
bases are related by 



Fa = \Y.^h Fl ° 



i -J 



F*3 



E(v 



1.1 



(49) 



A a are 3x3 Gell-Mann matrices [6]. The decomposition (]47p corresponds to 3x3 = 3 + 6 where 
the representation 3 is antisymmetric and corresponds to orbital angular momentum I 13 while 
the representation 6 is symmetric and correspond to the quadropole moment Q l] . In particular 
F 2 = mul\F 5 = -mu;l 2 ,F 7 = mcol 1 while F 1 ,F 3 ,F A ,F 6 and F 8 correspond to Q« e.g. F t = 



mux x + p < p = Qlf + Q~. Using the commutation relations, [I is 3 x 3 unit matrix] 



[x,P J 

[ Y'Y J 



il 



if, 



abc~ 



(50) 
(51) 



we obtain 



[F a ,F b ] =if abc muF c 



(52) 
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where a,b,c = 2, 5, 7 and this corresponds to 

[l\P]=ie^ k l k 



(53) 



On the other hand for a,b = 1,3,4,6,8 in Eq. (|5ip . c is restricted to 2,5,7 and /123 = 
1, (/i47, /246, /257, /34s) = 1/2, (/156, iW) = -1/2, (/458, /e7s) = \/3/2, all others are zero. 

Now for the relativistic harmonic oscillator where we have the Dirac Hamiltonian (3), with 
V(r) = ■^rnuP'r 2 , qaudrupole moment operator takes the form 

r« = (l + f3)f(r)Q^ + Q^ + 1 5 [(l-l3)g(r)Q^.p+(l + f3)^.pg(r)Q\ 
Since [ 7 5 S • p, Q P j ] = 0, 



ij 



(54) 



[IT, I*'] = 7 5 P-P, /(r)Q^-2(y + mMr)Q^]+7 5 /3[S-p, /(r)gi J -2(y + m) S (r)Q! 
+ (1 + /3)[F, Q^] + (l + /3)[p 2 , g(r)Q%] 

Thus [if, T*- 7 ] = gives [with Q l i = ?n 2 uj 2 (x l xi — |<5 u x 2 ) and similar expression for Qp] 



(55) 



/(r) = 2(V + m)g(r) 
1 

2?71 



s(r) 



with V(r) = \Mbj 2 r 2 . Thus 



r « = (i + /3 )-(y + m)Q^ + Q^ + — 7 5 (l-/3)£^£.p+(l + /3)£.pQ* 



2m 



In the Gell-Mann basis this becomes 



r o = (1 + $)-{V + m)FZ + F£ + -L 7 5 i 
m 2m 

Then using the commutation relations (50) and (51) 

[T a , r 6 ] = i/ a6c w ((1 + /3)(V + m)F c + ^ 7 5 1 1 



where F c = imcj[x T 4rp — p T ^fx\. Thus 
[T ,r 6 ] =ifab c u 



{H + m)F c + )p\l-p)[F c ,t-p\ 



(56) 



(57) 



(58) 



(59) 



(60) 



where a, 6 = 1, 3, 4, 6, 8 while c = 2, 5, 7, so that F c is essentially the orbital angular momentum I. 
Now using Eqs. (15) and (18) 



(H + m)L k = 7 5 S-p+(l + /3)(T/ + m) 



Z fc + -(l-/3)(£ fc 



^2~P 



(fl- + m )Z* + I 7 5 (l-/3)(Exp)* 



(J + ^ + ^l-^)- 

2 mw L 



Fc S • p 



(61) 
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where we have used ^ F c , £ • p ~ l k ,Y, ■ p 
where now F 2 = L 3 , F 5 = —L 2 , F 7 = L 1 and T a 



|(S x p) k , Thus from Eqs. (60) and (61) 



^Jabc-Fc 



a,b,c = 2, 5, 7 and 



Fa,Fb 



Fa,F b 



—, — =s =. Further 

yjmijj 2 (H +m) 
~ IJabc-Pc 



^Jabc'- c 



(62) 



(63) 



(64) 



where a = 2, 5, 7 and b, c = 1, 3, 4, 6, 8. 

The commutation relations (|621 I63l[64"|) generate the SU(3) algebra. Since S commutes with H 
as well as with all the above generators, the invariant algebra of the relativistic harmonic oscillator 
represented by the Dirac Hamiltonian (3) with V(r) = ^muj 2 r 2 is SU(3) <g> SU a (2). 

We now calculate the energy spectrum for which purpose we note from Eq. (58) that [a = 
1,3,4,6,8] 



T 2 = £r a r a = I (H 2 -m 2 ) 2 -mco 2 (H + m){L 2 + 3) 



(65) 



where we have used the definitions of F£ and F% given in Eq. (|45p . e.g. Ff = mu 2 x\X2 and 
Ff = p lP2 , Eqs. (11) and (19) with V = \muj 2 r 2 , 
Eq. (15). In terms of r a , Eq. (|65p takes the form 



¥f = pip 2 , Eqs. (11) and (19) with V = \mu?r 2 ', (1 + /3)(1 - 0) = 0, 75 (1 ± p) = (1 T PH and 



(H - m) 2 (H + m) - 9mw 2 



(66) 



3muj 2 J^r a r a + L 2 

But L 2 = J2b F b 2 , where 6 = 2,5, 7. Thus Eq. Q66) takes the form 

3?nw 2 f 2 = (H -m) 2 (H + m)-9mu 2 (67) 

where T 2 = Xla^a^a + Sb-T&rft is the invariant of the group SU(3) and as such is proportional to 
unit matrix. 

Hence Eq. (|67p gives the energy eigenvalues 

(E-m) 2 (E + m)-9muj 2 = Cmu} 2 (68) 

where C is to be fixed. This can be done in the following way. We take the non-relativistic limit 
of Eq. (18), H — > H non _ re i + m, which gives 



tj P . n T/ V . 2 2 

H non - r ei = h 2K = - Vmuj r 

2m 2m 



(69) 
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which, as is well known, gives the energy eigenvalues [note we have to replace to by \/2oj in the 
ordinary harmonic oscillator eigenvalues] 

e N = -^oo(2N + 3) (70) 



Now we take the non-relativistic limit of Eq. ()68[) . E — > £n + m, which gives 

moo 2 (2N + 3) 2 -9mu 2 = Cmco 2 (71) 

fixing C = AN(N + 3). Putting back in Eq. (|68p the energy eigenvalue equation becomes 

(E N -m) 2 (E N + m) = A(N + ^\ mu? (72) 



where N = 0, 1, • • • . This agrees with one obtained from the exact solutions of Dirac equation 



V. SUMMARY AND CONCLUSIONS 



Q. 



We have systematically reviewed the various dynamical symmetries of the Dirac Hamiltonian, 
clearly stating the conditions under which such symmetries hold. These symmetries include rel- 
ativistic spin (pseudo-spin) and orbital angular symmetries which hold when Dirac Hamiltonian 
with scalar Vs(r) and vector Vy{r) spherically symmetric potentials satisfy dVy/dr = ±dVs/dr or 
Vy = ±Vs + U. Here U is a constant potential but can be absorbed in redefinition of the mass. Then 
if V(r) = —a/r, as for the hydrogen atom or in QCD the Dirac Hamiltonian has 30(4)0 SU a (2) or 
equivalently SU M (2) <3 SU N (2) <3> SU a (2) symmetry. Here 577 CT (2) is the 577(2) group generated by 
the relativistic spin S l defined in Eq. (8) or Eq. (12). If V(r) = ^Mco 2 r 2 : , the harmonic oscillator 
potential, then the symmetry is SU(3) <S> SU a {2). We have used Dirac algebra of Dirac matrices 
and for the simple harmonic oscillation, Gell-Mann basis of 577(3), which is more transparent and 



a. 



simple (at least for physicist with particle physics background) compared to the basis used in 
We have also calculated energy spectrum in each case from group theoretical consideration, which 
agrees with the exact solution of Dirac equation in each case. We have also indicated physical 
systems where dynamical symmetries discussed above are possibly relevant. 
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